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ABSTRACT 

We investigate the evolution of high redshift seed black hole masses at late times 
and their observational signatures. The massive black hole seeds studied here form at 
extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating 
dark matter halos with seeds formed in this way, we follow the mass assembly of these 
black holes to the present time using a Monte-Carlo merger tree. Using this machinery 
we predict the black hole mass function at high redshifts and at the present time; the 
integrated mass density of black holes and the luminosity function of accreting black 
holes as a function of redshift. These predictions are made for a set of three seed 
models with varying black hole formation efficiency. Given the accuracy of current 
observational constraints, all 3 models can be adequately fit. Discrimination between 
the models appears predominantly at the low mass end of the present day black hole 
mass function which is not observationally well constrained. However, all our models 
predict that low surface brightness, bulgeless galaxies with large discs are least likely 
to be sites for the formation of massive seed black holes at high redshifts. The efficiency 
of seed formation at high redshifts has a direct influence on the black hole occupation 
fraction in galaxies at z = 0. This effect is more pronounced for low mass galaxies. 
This is the key discriminant between the models studied here and the Population III 
remnant seed model. We find that there exists a population of low mass galaxies that 
do not host nuclear black holes. Our prediction of the shape of the Mbh — o" relation 
at the low mass end is in agreement with the recent observational determination from 
the census of low mass galaxies in the Virgo cluster. 
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1 INTRODUCTION 

The demography of local galaxies suggests that most galax- 
ies host a quiescent supermassive black hole (SMBH) at 
the present time and the properties of the black hole 
are correlated with those of the host spheroid. In par- 
ticular, recent observational evidence points to the exis- 
tence of a tight correlation between the mass of the cen- 
tral blac k hole and the veloc i ty dispersion of the host 
spheroid iFTremaine et all |2002| : iFerrarese fc MerrittI l2000l : 
iGebhardt et al.l 2000l ) in nearby galaxies. This correlation 
strongly suggests coeval growth of the black hole and the 
stellar componen t via likely regulation of the gas supply 
in ga l actic nuclei llSilk fc ReeJl998l:lKauffmann fc Haehneltl 
I2OOOI : lKing||2003l : [Thompson et al.ll2005l ). 

Bla ck hole growth is b elieved to be powered by gas ac- 
cretion (|Lvnden-Belll Il969l ) and accreting black holes are 
detected as optically bright quasars. These optically bright 
quasars appear to exist out to the highest redshifts probed 



at the present time. Therefore, the mass build-up of SMBHs 
is likely to have commenced at extremely high redshifts [z > 
10). In fact, optically bright quasars have now been detected 
at z > 6 (e.g.. iFan et al..,2004 : Fan et al, 200^ in the Sloan 
Digital Sky Survey (SDSS). Hosts of hi gh redshift quasars 
are often strong sources of dust emission llOmont et al ]|200ll : 
ICox et aLlliooa : ICariUi et al.ll2002l : IWaher et al.ll2003l ) , sug- 
gesting that quasars were already in place in massive galax- 
ies at a time when galaxies were undergoing vigourous star 
formation. The growth spurts of SMBHs are also seen in 
the X-ray waveband. The integrated emission from these X- 
ray quasars generates the cosmic X-ray background (XRB) , 
and its spectrum suggests that most black-hole growth is 
obscured in optical wavelengths (iFabian fc Iwasawa 19991: 
Mushotzky et al.l I2OO0I: iHasinger et al.l I2001I ; iBarger et all 



20031 . I2OO5I : IWorslev et all bOOSl ). There exist examples of 



obscured black-hole growth in the form of 'Type-2' quasars, 
but their detected numbers are fewer than expected from 
models of the XRB. However, there is recent tantalizing ev- 
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idence from infra-r ed (IR) studies that dust-obscu red accre- 
tion is ubiquitous (jMartmez-Sansigre et al.ll2005l ). Current 
work suggests that while SMBHs might spend most of their 
hfetime in an optically dim phase, the bulk of mass growth 
occurs in the short-lived quasar stages. 

The assembly of BH mass in the Universe has been 
tracked using optical quasar activity. The current phe- 
nomenological approach to understanding the assembly of 
SMBHs involves optical data from both high and low red- 
shifts. These data are used as a starting point to con- 
struct a consistent picture that fits within the larger 



the Universe (|Haehnelt et al.l 


19981: iHaiman & Loebj 


1998 


Kauffmann & Haehneltl I2OO0I. 


2OO2I: IWvithe & Loebl 


2002 


Volonteri et all 20031: Di Matteo et all 2003f). 



Current modeling is grounded in the framework of the 
standard paradigm that involves the growth of structure 
via gravitational amplification of small perturbations in 
a CDM Universe — a model that has independent valida- 
tion, most recently from Wilkinson Microwave Anisotropy 
Probe (WMAP) measurements of the anisotropie s in the cos- 
mic m icrowave background (iSpergel et al. .200 j : iPage et al.l 
l2003l ). Structure formation is tracked in cosmic time by keep- 
ing a census of the number of collapsed dark matter halos 
of a given mass that form; these provide the sites for har- 
boring black holes. The computation of the mass function of 
dark matter halos i s don e using either the Press-Schechter 
l|Press fc Schechteil 19741) o r the extended Press-Schechter 
theory (|Lacev fc Col(j'l993h. or Montc-Carlo realizations of 
merger trees (K auffmaim & Haohnolt 2000; Volon teri et al.l 
l2003l : iBromiev et al.l 120041 ) or, in so me cases, directly from 
cosm ological N-body simulations l|Di Matteo et aL |2003| . 
I2OO5I ). 

In particular IVolonteri et al.l l|2003l ) have presented a 
detailed merger-tree based scenario to trace the growth of 
black holes from the earliest epochs to the present day. 
Monte-Carlo merger trees are created for present day halos 
and propagated back in time to a redshift of ~ 20. With the 
merging history thus determined, the initial halos at 2 ~ 20 
are then populated with seed black holes which are assumed 
to be remnants of the first stars that form in the Universe. 
The masses of these so-called Population III stars are not ac- 
curately known, however numerical simulatio ns by various 
groups (jAbel et al.ll2000l : iBromm et al.|[20o3 ) suggest that 
they are skewed to high masses of the order of a few hun- 
dred solar masses. Seeded with the end products of this first 
population, the merger sequence is followed and black holes 
are assumed to grow with every major merger episode. An 
accretion episode is assumed to occur as a consequence of 
every merger event. Following the growth and mass assem- 
bly of these black holes, it is required that the model is in 
consonance with the observed local A/bh — o" relation. The 
luminosity function of quasars is predicted by these models 
and can be compared to observations. Volonteri et al. find 
that not every halo at high redshift needs to be populated 
with a black hole seed in order to satisfy the observational 
constraints at 2: = 0. These models do not automatically re- 
produce the required abundance of supermassive black holes 
inferred to power the observed z > 6 SDSS quasars. In or- 
der to match the observation and produce SMBHs roughly 
1 Gyr after the Big Bang, it is required that black holes 
undergo brief, but extremely strong growth episodes dur- 



ing which the accre tion rate onto them is well in excess of 
the E ddington rate (jVolonteri fc Reesll2005l : [Begel man et al.l 
I2OO6I ). It is the existence of these SMBHs powering quasars 
at 2 > 6 that has prompted work on alternate channels to 
explain their mass build-up. 

In order to alleviate the problem of explaining the exis- 
tence of SMBHs in place by z ~ 6, roughly 1 Gyr after the 
Big Bang, in this paper we examine the possibility of using 
a well motivated high redshift seed black hole mass function 
as the initial black hole population at the highest redshifts. 
We investigate the effect of populating early dark matter 
halos with massive black hole se e ds predicte d in a model 
proposed bv lLodato fc NataraianI (|2006l . l2007l ). This model 
predicts a mass function for black holes that results from the 
direct collapse of pre-galactic gas discs. We study the impli- 
cations of the use of this seed mass function versus that of 
the Population HI remnants, in particular the difference in 
predictions at 2; = for the massive seed models. In Section 
2, we briefly outline the high redshift BH seed formation 
model, in Section 3 we evolve this model with redshift us- 
ing the merger-tree formalism developed by I Volonteri et al] 
(2003). The results are presented in Section 4, followed by a 
discussion of implications in the final section. 



2 BH SEED FORMATION MODEL 

In this paper, we track the formation of seed black holes 
in an ab-initio model and follow their mass assembly down 
to 2; = 0. This is done in two separate phases - starting 
with the high redshift seeds and tracing their subsequent 
growth. At high redshift (z > 15), we assume the inter- 
galactic medium has not been significantly enriched by met- 
als, and therefore the gas cooling timescale s are long. Under 
these conditions, many authors [Koushiappas et al. |2004| : 



iBegelman et al.ll2006l : lLodato fc Nataraianll2006l . l2007l ) have 
shown that pre-galactic discs can efflciently transport mat- 
ter into their innermost regions through the development 
and amplification of non-axisymmetric gravitational insta- 
bilities, often without fragmentation and star formation tak- 
ing place (see below). This is the main seed formation phase, 
wherein massive seeds with M « 10^ - 10^ Mq can form. At 
lower redshifts, cooling becomes more efficient and we as- 
sume that further accretion occurs via a 'merger driven sce- 
nario' described in more detail in Section 3. In this section we 
provide simple, analytical estimates of the amount of mass 
that we expect to be assembled in the form of massive BH 
seeds, based on the above scenario, as a function of the key 
dark matter h alo parameters. Here we ref er in particular to 
the model by iLodato fc NataraianI (|2006l ). who considered 
the evolution of pre-galactic discs by self-consistently tak- 
ing into account gravitational stability and fragmentation, 
thereby providing a detailed inventory of the fate of the gas. 

Consider a dark matter halo of mass M and virial tem- 
perature Tvir, containing gas mass Mgas = m^M (we also 
assume that the baryon fraction is roughly 5% implying 
nid — 0.05), of primordial composition, i.e. gas not enriched 
by metals, for which the cooling function is dominated by 
hydrogen. The other main parameter characterizing a dark 
matter halo that is relevant to the fate of the gas is its 



spin parameter A (= JhEj/^ /GM^' , where Jh is the total 
angular momentum and Eh is the binding energy) . The dis- 
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tribution of spin parameters for dark matter hales measured 
in numerical simulations is well fit by a lognormal distribu- 
tion in Aspin, with mean Aspin = 0.05 and standard deviation 
aA = 0.5: 



p(A)dA = 



exp 



ln2(A/A) 



2al 



dX 
T' 



(1) 



This function has been shown to provide a good fit 
to the N-body results of several i nvestigations (e.g. . 



Warren et al.ll 19921: ICole fc Lacevl[l99^ : iBuUock et allbOOll : 
van den Bosch et allbOoi T 

If the virial temperature of the halo Tvir > Tgas, the gas 
collapses and forms a rotationally supported disc. For low 
values of the spin parameter A the resulting disc can be com- 
pact and dense and is subject to gravitational instabilities. 
This occurs when the stability parameter Q defined below 
approaches unity: 



Q 



ttGE 



= ^/2. 



nGT.R' 



(2) 



where R is the cylindrical radial coordinate, E is the sur- 
face mass density, Cs is the sound speed, k = V^Vh/R is 
the epicyclic frequency, and Vh is the circular velocity of the 
disc (mostly determined by the dark matter gravitational 
potential). We have also assumed that at the relevant radii 
(« 10^ — 10^ pc) the rotation curve is well described by a 
flat Vh profile. We consider here the earliest generations of 
gas discs, which are of pristine composition with no met- 
als and therefore can cool only via hydrogen. In thermal 
equilibrium, if the formation of molecular hydrogen is sup- 
pressed, these discs are expected to be nearly isothermal 
at a temperat ure of a few thousand Ke lvin (here we take 
Tgas ~ 5000K. iLodato fc NataraianlbOO^ ). However, molec- 
ular hydrogen if present can cool these discs further down 
to temperatures of a few hundred Kelvin. The stability pa- 
rameter has a critical value Qc of the order of unity, below 
which the disc is unstable leading to the potential forma- 
tion of a seed black hole. The actual value of Qc essentially 
determines how stable the disc is, with lower Qc's imp lying 
more stable discs. It is well known since iToomrd l|l96A pro- 
posed this stability criterion, that for an infinitesimally thin 
disc to be stable to fragmentation, Qc = 1 for axisymmetric 
disturbances. The exact value of Qc under more realistic con- 
ditions is not well determined. Finite thickness efi'ects tend 
to stabilize the disc (reducing Qc), while on the other hand 
non-axisymmetric perturbations are in reality more unsta- 
ble (enhancing Qc) . Global, three-dimensiona l simulations 
of Keplerian discs (|Lodato fc Ricell2004l . [ioosi ) have shown 
that such discs settle down in a quasi-equilibrium configu- 
ration with Q remarkably close to unity, implying that the 
critical value Qc ~ 1. In this paper, we take Qc to be a free 
parameter and evaluate our results for a range of values. 

If the disc becomes unstable it develops non- 
axisymmetric spiral structures, which leads to an effective 
redistribution of angular momentum, thus feeding a growing 
seed black hole in the center. This process stops when the 
amount of mass transported to the center, A/bh, is enough 
to make the disc marginally stable. This can be computed 
easily from the stability criterion in eqn. ([2} and from the 
disc properties, determine d from the dark matter halo mass 
and angular momentum (|Mo et al.l Il998h . In this way we 



obtain that the mass accumulated in the center of the halo 
is given by: 



Mbh 



where 

Amax — 



rWdM 




8A 



mdQc \md J VTvir 



1/2 



: mdQc/8(md/jd)(Tvi,/Tga.) 



1/2 



A < An 
A > An 

(4) 



is the maximum halo spin parameter for which the disc is 
gravitationally unstable. Note that while Qc = 1 provides 
the benchmark for the onset of instability, non-axisymmetric 
and global instabilities can cause the disc to become unsta- 
ble for larger values of Qc. For this reason we investigate 
models with 1 < Qc ^ 3. 

The process described above provides a means to trans- 
port matter from a typical scale of a few hundred parsecs 
down to radii of a few AU. If the halo-disc system already 
possesses a massive black hole seed from a previous gener- 
ation, then this gas can provide a large fuel reservoir for 
its further growth. Note that the typical accretion rates im- 
plied by the above model are of the order of O.OIM© /yr, 
and are therefore sub-Eddington for seeds with masses of 
the order of 10^ Mq or so. If, on the other hand no black 
hole seed is present, then this large gas inflow can form a 
seed anew. The ultimate fate of the gas in this case at the 
smallest scales is more uncertain. One possibility, if the ac- 
cretion iute_Js_suffi£iently_large, has been described in de- 
tail bv lBegelman et al] (|2006f ). The infalling material likely 
forms a quasi-star, the core of which collapses and forms 
a BH, while the quasi-star keeps accreting and growing in 
mass at a rate which would be super-Eddington for the cen- 
tral BH. Alternatively, the gas might form a super-massive 
star, which would eventua lly collapse and form a black hole 
(|Shapiro fc Shibatal [2003 1 . There are no quantitative esti- 
mates of how much mass would ultimately end up collapsing 
in the hole. Thus, the black hole seed mass estimates based 
on eqn. |3] should be considered as upper limits. 

For large halo mass, the internal torques needed to re- 
distribute the excess baryonic mass become too large to be 
sustained by the disc, which then undergoes fragmentation. 
This occurs when the virial temperature exceeds a critical 
value Tniax, given by: 



Tgas 



\ md 1 + MBn/rridM 



2/3 



(5) 



where etc ~ 0.06 is a dimensionless parameter measuring the 
critical gravitatio nal torque above which the disc fragments 
(|Rice et al.ll2005h . 

To summarize, every dark matter halo is characterized 
by its mass M (or virial temperature Tvir) and by its spin 
parameter A. The gas has a temperature Tgas ~ 5000K. 
If A < Amax (see eqn. H} and Tvir < Tnax (eqn. O, then 
we assume that a seed BH of mass Mbh given by eqn. ([3} 
forms in the center. The remaining relevant parameters are 
rrid ~ 0.05, Oc = 0.06, and we consider three different values 
for Qc — 1.5, 2, 3, which will be referred to as (i) model A; 
(ii) model B and (iii) model C which correspond respectively 
to cases of increasing instability and therefore increasing ef- 
ficiency for the formation of seed black holes. What we in- 
vestigate in this paper are possible constraints/insights that 
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Figure 1. The probability of hosting a BH seed of any mass 
at ^ = 18 as a function of dark matter halo mass. The three 
curves refer to Qc = 1.5 (low efficiency case, solid line), Qc = 
2 (intermediate efficiency, short-dashed line) and Qc = 3 (high 
efficiency, long-dashed line). 

the measured mass function of supermassive black holes at 
z = can provide on the onset of instability and therefore on 
the efficiency of seed formation at extremely high redshifts. 

To give an idea of the efficiency of BH seed formation 
at high z within the present model, we plot in Fig. [1] as an 
example, the probability of forming a BH (of any mass) at 
z = 18, as a function of halo mass, for the three models. 
It can be seen that typically up to 10% of the haloes in the 
right mass range can form a central seed BH, the percentage 
rising to a maximum of ~25% for the high efficiency model 
C (highly unstable discs), and dropping to a maximum of 
~ 4% for the high stability and therefore low efficiency case 
(model A) . 



3 THE EVOLUTION OF SEED BLACK HOLES 

We follow the evolution of the MBH population result- 
ing from the seed formation process delineated above in a 
ACDM Universe. Our approach is similar to the one de- 
scribed in Volonteri, Haardt & Madau (2003). We simu- 
late the merger history of present-day halos with masses 
in the range lO" < M < 10^'' Mq starting from z = 20, 
via a Monte Carlo algorithm based on the extended Press- 
Schechter formalism. 

Every halo entering the merger tree is assigned a spin 
parameter according to eqn. [T] Recent work on the fate of 
halo spins during mergers in cosmological simulations has 
led to conflicting results: Vitvitska et al. (2002) suggest that 
the spin parameter of a halo increases after a major merger, 
and the angular momentum decreases after a long series 
of minor mergers; D'Onghia & Navarro (2007) find instead 
no significant correlation between spin and merger history. 
Given the unsettled nature of this matter, we adopt Occam's 



razor to guide us, and assume that the spin parameter of a 
halo is not modified by its merger history. 



When a halo enters the merger tree we assign seed 
MBHs by determining if the halo meets all the require- 
ments described in Section 2 for the formation of a cen- 
tral mass concentration. As we do not self-consistently trace 
the metal enrichment of the intergalactic medium, we con- 
sider here a sharp transition threshold, and assume that the 
MBH formation scenario suggested by Lod ato fc Nataraja n 
ceases at 2: ~ 15 (see also Sesana et al..2007i ; IVolonteril20CI7l ). 
At 2; > 15, therefore, whenever a new halo appears in the 
merger tree (because its mass is larger than the mass reso- 
lution), or a pre-existing halo modifies its mass by a merger, 
we evaluate if the gaseous component meets the conditions 
for efficient transport of angular momentum to create a large 
inflow of gas which can either form a MBH seed, or feed one 
if already present. 



The efficiency of MBH formation is strongly dependent 
on critical value of the Toomre parameter Qc, which sets 
the frequency of formation, and consequently the number 
density of MBH seeds. We investigate the inffuence of this 
parameter in the determination of the global evolution of the 
MBH population. Fig. [5] shows the number density of seeds 
formed in three models, with Qc = 1.5 (low efficiency model 
A), Qc = 2 (intermediate efficiency model B), and Qc = 3 
(high efficiency model C). The solid histograms show the 
total mass function of seeds formed by z = 15 when this 
formation channel ceases, while the dashed histograms refer 
to seeds formed at a specific redshift slice at z = 18. The 
number of seeds changes by about one order of magnitude 
from the least efficient to the most efficient model, consistent 
with the probabilities shown in Fig. 1. 



We assume that, after seed formation ceases, the 2 < 15 
population of MBHs evo lves according t o a "m erger driven 
scenario" , as described in lVolonteri et akl (|2006l ). We assume 
that during major mergers MBHs accrete gas mass that 
scales with the fifth power of the circular velocity ( or equiva- 
lently the velocity dispersion ac) of the host halo (|Ferraresel 
l2002f ). We thus set the final mass of the MBH at the end of 
the accretion episode to 90% of the mass predicted by the 
A/bh — correlation, assuming that the scaling does not 
evolve with redshift. Major mergers are defined as mergers 
between two dark matter halos with mass ratio between 1 
and 10. BH mergers contribute to the mass increase of the 
remaining 10%. 



In order to calculate the accreting black hole luminos- 
ity function and to follow the black hole mass growth during 
each accretion event, we also need to calculate the rate at 
which the mass, as estimated above, is accreted. This is as- 
sumed to scale with the Eddington rate for the MBH, and is 
based on the results of merger simulations , which heuristi- 
cally track accretion onto a central MBH (|Di Matteo et all 
l2005l : lHopk ins et al ]|2005ah . The time spent by a given simu- 
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Figure 2. Mass function of MBH seeds in three Q-models of tliat 
differ in seed formation efficiency. Left panel: Qc = 1.5 (the least 
efficient model A), middle panel: Qc = 2 (intermediate efficiency 
model B), right panel: Qc = 3 (highly efficient model C). Seeds 
form at 2 > 15 and this channel ceases a,t z = 15. The solid 
histograms show the total mass function of seeds formed hy z = 
15, while the dashed histograms refer to seeds formed at a specific 
redshift, z = 18. 



lated AGN at a given bolome tric luminosit; 
interval is approximated by l|Hopkins et al 



per lo garithmic 



20053) 



dt ( I r 



-1 f L 



(6) 



where tg ~ lO'' yr, and a = -0.95 + 0.32 log(Lpcak/lO^^L0). 
Here Lpcak is the luminosity of the AGN at the peak of its ac- 
tivity. Hopkins et al. (2006) show that approximating ipcak 
with the Eddington luminosity of the MBH at its final mass 
(i.e., when it sits on the Mbh — o^c relation) compared to 
computing the peak luminosity with eqn. (6) above gives 
the same result and in fact, the difference between these 2 
cases is negligible. Volonteri et al. (2006) derive the follow- 
ing simple differential equation to express the instantaneous 
accretion rate (/Edd,in units of the Eddington rate) for a 
MBH of mass Mbh in a galaxy with velocity dispersion ac'- 



dfEdd{t) /Edd (*) f eA^Eddc' 



dt 



lO^Lr. 



(7) 



where here t is the time elapsed from the beginning of the 
accretion event. Solving this equation gives us the instanta- 
neous Eddington ratio for a given MBH at a specific time, 

^ We convert accretion rate into luminosity assuming that the 
radiative efficiency equals the binding energy per unit mass of a 
particle in the last stable circular orbit. We associate the location 
of the last stable circular orbit to the spin of the MBHs, by self- 
consistently trac king the evolution of b lack hole spins throughout 
our calculations l| Volonteri et al.ll2005^ . We set 20% as the max- 
imum value of the radiative efficiency, corresponding t o a spin 
slight ly below the theoretical limit for thin disc accretion ijThornel 
1 1974) . 



and therefore we can self-consistently grow the MBH mass. 
We set the Eddington ratio /Edd ~ at t = 0. This same 
type of accretion is assumed to occur, at z > 15, following a 
major merger in which a MBH is not fed by disc instabilities. 

In a hierarchical Universe, where galaxies grow by 
mergers, MBH mergers are a natural consequence, and 
we trace their contribution to the evolving MBH popula- 
tion (cfr. Sesana et al. 2007 for details on the dynami- 
cal modeling). During the final phases of a MBH merger, 
emission of gravitational radiation drives the orbital de- 
cay of the binary. Recent numerical relativity simulations 
suggest that merging MBH binaries might be subject to 
a large "gravitational recoil": a genera l-relativistic effect 
l|Fitchettl il983l : iRedmount fc Ree j [l989l ') due to the non- 
zero net linear momentum carried away by gravitational 
waves in the coalescence of two unequal mass black holes. 
Radiation recoil is a strong field effect that depends on 
the lack of symmetry in the system. For merging MBHs 
with high spin, in particular orbital configurations, the re- 
coil velocit y can be as high as a few th ousands of kilometers 
per second ('Campanolli ct al."2007a''b'; iGonzalez et al.ir2007l : 
Herrmann c^al. 2007; Schnittman 20071)" Here, we aim to 
determine the characteristic features of the MBH popula- 
tion deriving from a specific seed scenario, and its signature 
in present-day galaxies, we study the case without gravita- 
tional recoil. We discuss this issue further in section 4. 



4 RESULTS 

Detection of gravitati onal waves from se eds merging at the 
redshift of formation (jSesana et al.|[2007l ) is probably one of 
the best ways to discriminate among formation mechanisms. 
On the other hand, the imprint of different formation sce- 
narios can also be sought in observations at lower redshifts. 
The various seed formation scenarios have distinct conse- 
quences for the properties of the MBH population a,t z — 0. 
Below, we present theoretical predictions of the various seed 
models for the properties of the local SMBH population. 



4.1 Superrticissive black holes in dwarf galaxies 

The repercussions of different initial efficiencies for seed for- 
mation for the overall evolution of the MBH population 
stretch from high-redshift to the local Universe. Obviously, 
a higher density of MBH seeds implies a more numerous 
population of MBHs at later times, which can produce ob- 
servational signatures in statistical samples. More subtly, 
the formation of seeds in a ACDM scenario follows the cos- 
mological bias. As a consequence, the progenitors of massive 
galaxies (or clusters of galaxies) have a higher probability of 
hosting MBH seeds (cfr. Madau & Rees 2001). In the case 
of low-bias systems, such as isolated dwarf galaxies, very 
few of the high-2 progenitors have the deep potential wells 
needed for gas retention and cooling, a prerequisite for MBH 
formation. We can read off directly from Fig. [T] the average 
number of massive progenitors required for a present day 
galaxy to host a MBH. In model A, a galaxy needs of order 
25 massive progenitors (mass above ~ 10^ M©) to ensure a 
high probability of seeding within the merger tree. In model 
C, instead, the requirement drops to 4 massive progenitors. 
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Figure 3. The Afbh— velocity dispersion (cc) relation at 2: = 0. Every circle represents the central MBH in a halo of given (Tc. 
Observational data are marked by their quoted errorbars, both in CTc, and in M^h (Tremaine et al. 2002). Left to right panels: Qc = 1.5, 
Qc = 2, Qc = 3, Population III star seeds. Top panels: fraction of galaxies at a given velocity dispersion which do not host a central 
MBH. 



increasing the probability of MBH formation in lower bias 
halos. 

The signature of the efficiency of the formation of MBH 
seeds will consequently be stronger in isolated dwarf galax- 
ies. Fig. [3] (bottom panel) shows a comparison between the 
observed Mbh — cr relation and the one predicted by our 
models (shown with circles), and in p articular, from left to 
right, the th ree models based on the iLodato fc NataraianI 
l|2006l . l2007l 'l seed masses with Qc = 1.5, 2 and 3, and a 
fourth model based on lower-mass Population HI star seeds. 
The upper panel of Fig. [3] shows the fraction of galaxies 
that do not host any massive black holes for different ve- 
locity dispersion bins. This shows that the fraction of galax- 
ies without a MBH increases with decreasing halo masses 
at z = 0. A larger fraction of low mass halos are devoid of 
central black holes for lower seed formation efficiencies. Note 
that this is one of the key discriminants between our models 
and those seeded with Population HI remnants. As shown 
in Fig. 3, there are practically no galaxies without central 
BHs for the Population III seeds. 

It is interesting to note that our model predictions are 
in very good agreement with the recent HST ACS cen- 
sus of black holes i n low mass galaxies in the Vi rgo clus- 
ter. [^rrareseeEai] l|2006l '): IWehner fc Harris (120061 ) suggest 
that below a transition galaxy mass (~ 10^° Mq) a central 



massive black hole seems to be replaced by a nuclear star 
cluster. Although no definite proof that Virgo dwarfs are in- 
deed MBH-less, the above results imply that MBHs are more 
common in large galactic systems. Our models also indicate 
that a minimum velocity dispersion exists, below which the 
probability of finding a central object is very low. 

We make quantitative predictions for the local occupa- 
tion fraction of MBHs. Our model A predicts that below 
(jc ~ 60 kms~^ the probability of a galaxy hosting a MBH is 
negligible. With increasing MBH formation efficiencies, the 
minimum mass for a galaxy that hosts a MBH decreases, 
and it drops below our simulation limits for model C. On 
the other hand, models based on lower mass Population III 
star remnant seeds, predict that massive black holes might 
be present even in low mass galaxies. 

We note here that in our investigation we have not in- 
cluded any mechanism that could further lower the occupa- 
tion fraction of MBHs (e.g., gravitational recoil, three-body 
MBH interactions). For any value of Qc the occupation frac- 
tion computed above is therefore an upper limit. 

Although there are degeneracies in our modeling (e.g., 
between the minimum redshift for BH formation and insta- 
bility criterion), the BH occupation fraction, and the masses 
of the BHs in dwarf galaxies are the key diagnostics. In local 
observations, the clearest signatures of massive seeds com- 



Evolution of massive black hole seeds 7 



-2 




a 10 12 14 8 10 12 14 
log(L/Lo) 



Figure 4. Bolometric luminosity functions at different redshifts. 
All 3 models match the observed bright end of the LF at high red- 
shifts and predict a steep slope at the faint end down to 2; = 1. 
The 3 models are not really distinguishable with the LF. How- 
ever at low redshifts, for instance at z = 0.5, all 3 models are 
significantly flatter at both high and low luminosities and do not 
adequately match the current data. As discussed in the text, the 
LF is strongly determined by the accretion prescription and what 
we see here is simply a reflection of that fact. 



pared to Population III remnants, would be a lower limit of 
order the typical mass of seeds (Fig. [2]) to the mass of MBHs 
in galaxy centers, as shown in Fig. [S] An additional caveat 
worth mentioning is the possibility that a galaxy is devoid 
of a central MBH because of dynamical ejections (due to ei- 
ther the gravitational recoil or three-body scattering). The 
signatures of such dynamical interactions should be more 
prominent in dwarf galaxies, but ejected MBHs would leave 
observational signatures on their hos t s (Gi ilteki n et al. in 
prep.) . On top of that, ISchnittmanI (|2007t ) and IVolonteril 
12003) agree in considering the recoil a minor correction to 
the overall distribution of the MBH population at low red- 
shift (cfr. figure 4 in Volonteri 2007). 

Additionally, as MBH seed formation requires halos 
with low angular momentum (low spin parameter), we en- 
visage that low surface brightness, bulge-less galaxies with 
high spin parameters (i.e. large discs) are systems where 
MBH seed formation is less probable. Furthermore, bul- 
geless galaxies are believed to preferentially have quieter 
merger histories and are unlikely to have experienced any 
major merger, which could have brought in a MBH from a 
companion galaxy. The (possible) absence of a MBH in M33 
hence arises naturally (Merritt, Ferrarese & Joseph 2001; 
Gebhardt et al. 2001) in our model. 



^ This prediction is pertinent to all models relying on gravita- 
tional instabilities triggered in low spin parameter halos. 



4.2 The luminosity function of accreting black 
holes 

Turning to the global properties of the MBH population, as 
suggested by Yu & Tremaine (2002) the mass growth of the 
MBH population at 2 < 3 is dominated by the mass accreted 
during the bright epoch of quasars, thus washing out most 
of the imprint of initial conditions. This is evident when we 
compute the luminosity function of AGN. Clearly the de- 
tailed shape of the predicted luminosity function depends 
most strongly on the accretion prescription used. With our 
assumption that the gas mass accreted during each merger 
episode is proportional to , we find that distinguishing be- 
tween the various seed models is difficult. As shown in Fig. U 
all 3 models reproduce the bright end of the observed bolo- 
metric LF (jHopkins et a l. 2007) at higher redshifts (marked 
as the solid curve in all the panels) , and predict a fairly steep 
faint end that is as yet undetected. All models fare less well 
at low redshift, shown in particular at z — 0.5. This could be 
due to the fact that we have used a single accretion prescrip- 
tion to model growth through epochs. On the other hand, 
the decline in the available gas budget at low redshifts (since 
the bulk of the gas has been consumed by this epoch by star 
formation activity) likely changes the radiative efficiency of 
these systems. Besides, observations suggest a sharp decline 
in the number of actively accreting black holes at low red- 
shifts across wave-lengths, produced most probably due to 
changes in the accretion flow as a result of change in the ge- 
ometry of the nuclear regions of galaxies. In fact, all 3 of our 
models under-predict the slope at the faint end. There are 
three other effects that could cause this flattening of the LF 
at the faint end at low redshift for our models: (i) not having 
taken into account the fate of on-going merging and the fate 
of satellite galaxies; (ii) the number of realizations generated 
and tracked is insufBcient for statistics, as evidenced by the 
systematically larger errorbars and (iii) more importantly, 
it is unclear if merger-driven accretion is indeed the trigger 
of BH fueling in the low redshift universe. We note that the 
3 massive seed models and Population III seed model can- 
not be discriminated by the LF at high redshifts. Models B 
and C are also in agreement viz-a-viz the predicted BH mass 
function at z = 6 (see Fig. 2), even assuming a very high 
radiative efficiency (up to 20%), while model A might need 
less severe assumptions, in particular for BH masses larger 
than lO'^ M0. 



4.3 Comoving mass density of black holes 

Since during the quasar epoch MBHs increase their mass by 
a large factor, signatures of the seed formation mechanisms 
are likely more evident at earlier epochs. We compare in 
Fig. [5] the integrated comoving mass density in MBHs to 
the expectations from Soltan-type arguments (F. Haardt, 
private communication) , assuming that qua sars are powered 
by radiatively efficient flows (for details, seelYu fc Tremainel 
I2OO2I : lElvis et al.ll2002l : [Marconi et all 120041 ) . While during 
and after the quasar epoch the mass densities in models A, 
B, and C differ by less than a factor of 2, at z > 3 the 
differences become more pronounced. 

A very efficient seed MBH formation scenario can lead 
to a very large BH density at high redshifts. For instance, 
the highest efficiency model C with Qc = 3, the integrated 
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Figure 5. Integrated black hole mass density as a function of 
redshift. Solid lines: total mass density locked into nuclear black 
holes. Dashed lines: integrated mass density accreted by black 
holes. Models based on BH remnants of Population III stars (low- 
est curve), Qc = 1.5 (middle lower curve), Qc = 2 (middle up- 
per curve), Qc = 3 (upper curve). Shaded area: constraints from 
Soltan-type arguments, where we have varied the radiative ef- 
ficiency from a lower limit of 6% (applicable to Schwarzschild 
MBHs, upper envelope of the shaded area), to about 20% (Wang 
et al. 2006). All 3 massive seed formation models are in comfort- 
able agreement with the mass density obtained from integrating 
the optical luminosity functions of quasars. 

MBH density at 2: = 10 is already ~ 25% of the density 
at 2; = 0. The plateau at 2; > 6 is due to our choice of 
scaling the accreted mass with the 2 = Mbh— velocity dis- 
persion relation. Since in our models we let MBHs accrete a 
mass which scales with the fifth power of the circular veloc- 
ity of the halo, the accreted mass is a small fraction of the 
MBH mass (see the discussion in lMaruUi et al.ll2006l ). and 
the overall growth remains small, as long as the mass of the 
seed is larger than the accreted mass, which, in our assumed 
scaling, happens when the mass of the halo is below a few 
times 10^" M0. The comoving mass density, an integral con- 
straint, is reasonably well determined out to 2: = 3 but is 
poorly known at higher redshifts. All models appear to be 
satisfactory and consistent with current observational limits 
(shown as the shaded area). 

4.4 Black hole mass function at 2 = 

One of the key diagnostics is the comparison of the mea- 
sured and predicted BH mass function at 2 = for our 3 
models. In Fig. |6l we show (from left to right, respectively) 
the mass function predicted by models A, B, C and Pop- 
ulation III remnant seeds compared to that obtained from 
measurements. The histograms show the mass function ob- 
tained with our models (where the upper histogram includes 
all the black holes while the lower one only includes black 
holes found in central galaxies of halos in the merger-tree 
approach). The two lines are two different estimates of the 



observed black hole mass function. In the upper one, the 
measured velocity dispersion function for nearby late and 
early-typo galaxies from the SDSS survey l|Bernardi et al.l 
[2003. ; Sheth et al. .2003 ') has been convolved with the mea- 
sured Mbh — o" relation. We note here that the scatter in the 
Afbh — cr relation is not explicitly included in this treatment, 
however the inclusion of the scatter is likely to preferen- 
tially affect the high mass end of the BHMF, which pro- 
vides stronger constraints on the accretion hi stories rather 
than the seed masses . It has been argued (e.g ., iTundo et al.l 
I2OO7I : iBernardi et ahllioOTl : iLauer et allbOOTi ) that the BH 
mass function differs if the bulge mass is used instead of 
the velocity dispersion in relating the BH mass to the host 
galaxy. Since our models do not trace the formation and 
growth of stellar bulges in detail, we are restricted to using 
the velocity dispersion in our analysis. 

The lower dashed curve is an alternate theoretical es- 
timate of the BH mass f unction derived usin g the Press- 
Schechter formalism from IJenkins et all l|200ll ) m conjunc- 
tion with the observed Mbh — cr relation. Selecting only the 
central galaxies of halos in the merger-tree approach adopted 
here (lower histograms) is shown to be fairly equivalent to 
this analytical estimate, and this is clearly borne out as is ev- 
ident from the plot. When we include black holes in satellite 
galaxies (upper histograms, cfr. the discussion in Volonteri, 
Haardt & Madau 2003) the predicted mass function moves 
towards the estimate based on SDSS galaxies. The higher 
efficiency models clearly produce more BHs. At higher red- 
shifts, for instance at 2 = 6, the mass functions of active 
MBHs predicted by all models are in very good agreements 
in particular for BH masses larger than 10® M©, as it is 
the growth by accretion that dominates the evolution of the 
population. At the highest mass end (> 10** M0) model A 
lags behind models B and C, although we stress once again 
that our assumptions for the accretion process are very con- 
servative. 

The relative differences between models A, B, and C at 
the low-mass end of the mass function, however, are gen- 
uinely related to the MBH seeding mechanism (see also 
Figs. |3] and In model A, simply, fewer galaxies host 
a MBH, hence reducing the overall number density of black 
holes. Although our simplified treatment does not allow ro- 
bust quantitative predictions, the presence of a "bump" at 
2 = in the MBH mass function at the characteristic mass 
that marks the peak of the seed mass function (cfr. Fig. [2]) 
is a sign of highly efficient formation of massive seeds (i.e., 
much larger mass with respect, for instance. Population 111 
remnants). The higher the efficiency of seed formation, the 
more pronounced is the bump (note that the bump is most 
prominent for model C). Since current measurements of 
MBH masses extend barely down to Mbh ~ 10® M© , this 
feature cannot be observationally tested with present data 
but future campaigns, with the Giant Magellan Telescope or 
JWST, are likely to extend the mass function measurements 
to much lower black hole masses. 



5 DISCUSSION AND IMPLICATIONS 

In this paper, we have investigated the role that the choice 
of the initial seed black hole mass function at high redshift 
(2 ~ 18) plays in the determination of observed properties of 



Evolution of massive black hole seeds 9 




456789456789456789456789 
lg(MB„(Mj) lg(M3„(Mj) lg(M3„(Mj) lg(M,„(Mg)) 



Figure 6. Mass function of blaclt holes at z=0. Histograms rep- 
resent the results of our models, including central galaxies only 
(lower histograms with errorbars), or including satellites in groups 
and clusters (upper histograms). Left panel: Qc = 1.5, mid-left 
panel: Qc = 2, mid-right panel: Qc = 3, right panel: models based 
on BH remnants of Population III stars. Upper dashed line: mass 
function derived from combining the velocity dispersion function 
of Sloan galaxies (Sheth et al. 2003, where we have included the 
late- type galaxies extrapolation), and BH mass- velocity disper- 
sion correlation (e.g., Tremaine et al. 2002). Lower dashed line: 
mass function derived using the Press-Schechter formalism from 
Jenkins et al. (2001) in conjunction with the Mbh ~ relation 
(Ferrarese 2002). 



local quiescent SMBHs. While the errors on mass determi- 
nations of local black holes are large at the present time, 
definite trends with host galaxy properties are observed. 
The tighest correlation appears to be between the BH mass 
and the velocity dispersion of the host spheroid. Starting 
with the ab-initio black hole seed mass function computed 
in the context of direct formation of central objects from the 
collapse of pre-galactic discs in high redshift halos, we fol- 
low the assembly history to late times using a Monte-Carlo 
merger tree approach. Key to our calculation of the evolu- 
tion and build-up of mass is the prescription that we adopt 
for determining the precise mass gain during a merger. Mo- 
tivated by the phenomenological observation of AIbh oc V^, 
we assume that this proportionality carries over to the gas 
mass accreted in each step. With these prescriptions, a range 
of predictions can be made for the mass function of black 
holes at high and low z, and the integrated mass density of 
black holes, all of which are observationally determined. We 
evolve 3 models, designated model A, B and C which corre- 
spond to increasing efficiencies respectively for the formation 
of seeds at high redshift. These models are compared to one 
in which the seeds are remnants of Population III stars. 

It is important to note here that one major uncertainty 
prevents us from making more concrete predictions: the un- 
known metal enrichment history of the Universe. Key to the 
implementation of our models is the choice of redshift at 
which massive seed formation is quenched. The direct seed 



formation channel described here ceases to operate once the 
Universe has been enriched by metals that have been synthe- 
sized after the first generation of stars have gone supernova. 
Once metals are available in the Inter-Galactic Medium gas 
cooling is much more efficient and hydrogen in either atomic 
or molecular form is no longer the key player. In this work 
we have assumed this transition redshift to be 2; = 15. The 
efficiency of MBH formation and the transition redshift are 
somehow degenerate (e.g., a model with Q = 1.5 and en- 
richment redshift 2: = 12 is halfway between model A and 
model B); if other constraints on this redshift were available 
we could considerably tighten our predictions. 

Below we list our predictions and compare how they 
fare with respect to current observations. The models in- 
vestigated here clearly differ in predictions at the low mass 
end of the black hole mass function. With future observa- 
tional sensitivity in this domain, these models can be dis- 
tinguished. 

Our model for the formation of relatively high-mass 
black hole seeds in high-z halos has direct influence on the 
black hole occupation fraction in galaxies at z = 0. This 
effect is more pronounced for low mass galaxies. We find 
that a significant fraction of low-mass galeixies might not 
host a nuclear black hole. This is in very good agreement 
with the shape of the Mbh — cr relation determined recently 
from an observational census (an HST ACS survey) of low 
mass galaxies in the Virgo cluster reported by Ferrarese et 
al. (2006). 

The models studied here (with different black hole seed 
formation efficiency) are distinguishable at the low mass end 
of the BH mass function, while at the high mass end the 
effect of initial seeds appears to be sub-domi nant. While 
current data in the low m ass regime is scant (iBarth et al.l 
I2OO4I : ICreene fc HollioOTi ). future instruments and surveys 
are likely to probe this region of parameter space with sig- 
nificantly higher sensitivity. 

All our models predict that low surface brightness, 
bulge-less galaxies with high spin parameters (i.e. large 
discs) are systems where MBH formation is least probable. 

One of the key caveats of our picture is that it is unclear 
if the differences produced by different seed models on ob- 
servables at 2: = might be compensated or masked by BH 
fueling modes at earlier epochs. There could be other chan- 
nels for BH growth that dominate at low redshifts like mi- 
nor mergers, dynamical instabilities, accretion of molecular 
clouds, tidal disruption of stars. The decreased importance 
of the merger driven scenario is patent from observations of 
low-redshift AGN, which are f or the large majori ty hosted 
by undisturbed galaxies (e.g., IPierce et al.l 120071 . and ref- 
erenc es therein) in low-density environments (e.g., iLi et al.l 
I2OO6I ). However, the feasibility and efficiency of some alter- 
native channels are still to be proven (for example, about 
the efficiency of feeding from large scale instabilities see 
discussion in iKing fc Pringld l2007i: IShlosman et all Il989l : 



lGoodmanll2003l : ICollin fc ZahnlllQQgf T In any event, while 
these additional channels for BH growth can modify the de- 
tailed shape of the mass function of MBHs, or of the lumi- 
nosity function of quasars, they will not create a new MBH. 
The occupation fraction of MBHs (see figure 3) is therefore 
largely independent of the accretion mechanism and a true 
signature of the formation process. 

To date, most theoretical models for the evolution of 
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MBHs in galaxies do not include how MBHs form. This work 
is a first analysis of the observational signatures of massive 
black hole formation mechanisms in the low redshift uni- 
verse, complementary to the investigation by Sesana et al. 
(2007), where the focus was on detection of seeds at the 
very early times where they form, via gravitational waves 
emitted during MBH mergers. We focus here on possible 
dynamical signatures that forming massive black hole seeds 
carries over to the local Universe. We believe that the signa- 
tures of seed formation mechanisms will be far more clear if 
considered jointly with the evolution of the spheroids that 
they host. The mass, and especially the frequency, of the 
forming MBH seeds is a necessary input when investigat- 
ing how the feedback from accretion onto MBHs influences 
the host galaxy, and is generally introduced in numerical 
models using extremely simplified, ad hoc prescriptions (e .g., 
Springcl ct al. 2005; Di Mattco et al. 2005; Hopkins o t al. 
2006i : ICroton et al.ll2005l : ICattaneo et al.ll2006l : IBower et al. 
2009 ). Adopting more detailed models for black hole seed 
formation, as outlined here, can in principle strongl y af- 
fect such results. For instance. iKauffmann et ahl l|2004) find 
that AGN activity is typically confined to galaxies with 
M > 10^" M0 . If we consider the occupation fraction of 
MBHs in such galaxies, we find that it differs by a large fac- 
tor between models A and C, being of order 10% in the low 
efficiency model (at 2; ~ 1 — 4) and 50% or higher in model 
C. Consequently, the possibility of AGN feedback and its 
effect on the host would be selective in the former case, or 
widespread in the latter case. Adopting sensible assump- 
tions for the masses, and frequency of MBH seeds in models 
of galaxy formation is necessary if we want to understand 
the symbiotic growth of MBHs and their hosts. 
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